Enhancement of critical current density and vortex activation energy in 
proton-irradiated Co-doped BaFe 2 As 2 
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The effect of proton irradiation in Ba(Feo.93Coo.o7)2As2 single crystals is reported. We analyze 
temperature dependence of current density and normalized flux relaxation rate in the framework of 
collective creep model. Glassy exponent and barrier height for flux creep are directly determined by 
Maley's method. Our model functions for barrier height and critical current density in the absence of 
flux creep are explained by the superposition of 8T C - and <5/-pinning. We also approach true critical 
current density by means of generalized inversion scheme, and the obtained result is in reasonable 
agreement with our model function. Proton irradiation effect on temperature dependence of current 
density and normalized relaxation rate can be summarized as doubling of barrier height at the 
beginning of flux creep. 

PACS numbers: 74.25.Wx, 74.25. Uv, 74.25. Sv, 74.70.Xa 



In high-temperature superconductors, many interest- 
ing phenomena in vortex dynamics are discovered such 
as giant-flux creep, thermally activated flux flow, and 
the theories to describe them have been elaborated in 
the last decades i Especially, collective pinning with 
weak pinning potential by the quenched disorder and 
collective creep of vortex bundles give rise to intrigu- 
ing experimental results, such as "plateau" observed 
in temperature dependent normalized relaxation rate 
(S = |dlnM/dlnt|)j2 in contrast to linear increase with 
temperature predicted by Anderson-Kim model in low- 
temperature superconductors. Recently discovered iron- 
based superconductors (IBSs) have relatively high critical 
temperature (T c ) and large critical current density (J c ). 
Besides, magnetization hysteresis loop in this system is 
quite similar to that in Y-Ba-Cu-0, and magnetic relax- 
ation measurements have revealed that IBSs also show 
giant-flux creep, which implies that IBSs and cuprate su- 
perconductors share common vortex physics. Moreover, 
how introduction of artificial pinning center affects flux 
dynamics or J c is also interesting.—' 4 In Y-Ba-Cu-0, J c 
is enhanced and glassy behavior remains basically the 
same after proton (H + ) irradiation, which is known to 
introduce point defects. This is also expected in IBSs. 
In fact, Haberkorn et al. recently reported that H + - 
irradiation does not affect to H — T/T c phase diagram in 
Ba(Feo. 925C00. 075)2 As2-— It is important to clarify how 
vortex dynamics is affected by H + -irradiation in IBSs. 

As we mentioned above, IBSs are suitable candidates 
to check whether glassy behavior of vortices is univer- 
sal in all high-temperature superconductors. Since this 
system is twin free and less anisotropic, it enables us 
to discuss intrinsic pinning and dynamic properties of 
vortices without complication. However, it is difficult to 
synthesize large and clean crystals especially in LnFeAsO 



(Ln is lanthanoid, so-called Till' system).— This pre- 
vents us from discussing vortex dynamics due to strong 
inhomogeneities unless we use local probes. Since high- 
quality single crystals are readily available in so-called 
'122' crystals AEFe 2 As 2 (AE : alkaline earths)^ it is 
possible to discuss details of vortex pinning also in IBSs 
with global magnetic measurements. Actually, homoge- 
neous flow of superconducting current in this system has 
been confirmed by magneto-optical measurement^ This 
is why we choose optimally Co-doped BaFe2As2 single 
crystal. 

In this paper, we report the effect of proton (H + ) ir- 
radiation in Ba(Feo.93Coo.o7)2As2 single crystals. We 
analyze it in the framework of collective creep theory 
with temperature dependent shielding current J and S. 
Glassy exponent and barrier height for flux creep are di- 
rectly determined. Our model functions for barrier height 
and critical current density in the absence of flux creep 
are explained by the superposition of 5T C - and JZ-pinning. 
We also approach the true critical current density by 
means of generalized inversion scheme, and the obtained 
result is in reasonable agreement with our model func- 
tion. Effects of H + irradiation on J(T) and S(T) can be 
summarized as doubling of barrier height at the begin- 
ning of flux creep. 

Optimally Co-doped BaFe2As2 single crystals were 
grown by FeAs/CoAs self- flux method. Fundamen- 
tal properties of this system have been reported else- 
wnerej 9 i io samples are cleaved to be thin plates 
with thickness less than ~ 20 /im. This value is much 
smaller than the projected range of 3 MeV H + for 
Ba(Feo.93Coo.o7)2 AS2 of ~ 50 //m, calculated by the stop- 
ping and range of ions in matter-2008»i 3 MeV H + ir- 
radiation was performed parallel to c-axis at 40 K at 
NIRS-HIMAC. The total dose of the measured sample is 
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1.2 x 10 16 cm" 2 . 1 2 Magnetization is measured in a com- 
mercial SQUID magnetometer (MPMS-XL5, Quantum 
Design) with a magnetic field parallel to c-axis. Main 
features of vortex dynamics in H + -irradiated sample an- 
alyzed in this work have been briefly reported in Rcf. 
Ila . To clarify the effect of H + irradiation on the vortex 
system, we also measured a pristine sample similar to 
the "unirradiated" sample in Ref. Current densities 
calculated by the Bean model are denoted as J s in field- 
sweep measurements (s is an abbreviation for sweep), or 
simply J in relaxation measurements. 

Figure [T] shows magnetic field dependence of J s in (a) 
pristine and (b) H + -irradiated Ba(Fe .93Co .o7)2As2. It 
is obvious that H + irradiation enhances J s from 1 x 10 6 
A/cm 2 to 2.5 x 10 6 A/cm 2 at 2 K under zero-field. In 
pristine sample, J s is nearly constant below 1 kOe, fol- 
lowed by power law decay H~ a at a field of 2-10 kOe with 
a ~ 0.5. As discussed by van der Beek et al, these be- 
haviors at low-fields are attributed to sparse strong-point 
pinning centers, as in the case of Y-Ba-Cu-0 filmsj 15 i 16 
At a glance, it seems inappropriate to analyze it in the 
framework of collective-pinning-collective-creep and vor- 
tex glass theory. However, since strong-point-pinning 
contribution for temperature (T) dependence of current 
density (J) is smaller than weak-collective pinning (see 
Fig. 6 or 9 in Ref. 1151 ), we can approximate J(T) only 
by the contribution from collective creep (pinning). In 
H + -irradiated sample, it is basically the same as the 
pristine one, although there is a wide-crossover region 
with a ~ 0.3 between low-field plateau and H~ 5 region. 
Such a weak field dependence has been also observed in 
YBa2Cu3C>7 and YBa2Cu40s films by Griessen et at, 
and they concluded that single vortex creep is achieved 
in this region*^ 

To elucidate the vortex dynamics, it is important 
to measure (static) magnetic relaxation rate 5* = 
|dlnA//dlni| in both samples, where M is magnetization, 
t is time from the moment when the critical state is pre- 
pared. In order to discuss temperature dependence of 
vortex dynamics, we have to fix magnetic field. How- 
ever, as we mentioned above, there is a strong pinning 
background, so that we have to select a field where field 
dependence of J s is similar for all temperatures to exclude 
field-dependent strong-point-pinning effect. Besides, we 
should carefully keep away from fish-tail effect with non- 
monotonic field dependence of J s at high fields and self- 
field effect at low fields, which disturb direct extraction of 
typical parameters for vortex dynamics. Based on these 
consideration, we select H = 5 kOe in both samples, 
shown as vertical broken lines in Fig. [T] Insets of Fig. [5] 
show temperature dependence of S. According to collec- 
tive pinning theory, — this is described as 
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U + LiT\n(t/t eS )' 



(1) 




where Uq is temperature-dependent flux activation en- 
ergy in the absence of flux creep, /i > is a glassy 
exponent for elastic creep, and i e ff is effective hopping 



/y(kOe) 

FIG. 1. Field dependence of J s in (a)pristine and (b)H + irra- 
diated Ba(Feo.93Coo.o7)2As2 at several temperatures. Dotted 
line on • (2 K) shows power-law decay of (a) H~ 5 and (b) 
H~ 3 , respectively. Vertical line indicates the field where we 
discuss the vortex dynamics, H = 5 kOe. 



attempt time. One of the most remarkable results ex- 
tracted from this equation is the prediction of plateau 
in the intermediate temperature range if Uq <C T. 2 The 
value of plateau S ~ l/(/xln(i/i ff)) falls in the range of 
0.02-0.04 theoretically, which has been confirmed in Y- 
Ba-Cu-Oi^ The inset of Fig. [Ha) is consistent with this 
behavior quantitatively, as observed in other IBSs. 1 — i 19 i 20 
This proves the validity of applying collective pinning 
theory to IBSs. Here we emphasize that it is quite im- 
portant to determine the value of fi in discussing vor- 
tex dynamics, since [i includes information on the size 
of vortex bundle in collective pinning theory. In three- 
dimensional system, it is predicted as fi — 1/7,3/2,7/9 
for single- vortex, small-bundle, and large-bundle regime, 
respectively^ Inverse current-density dependence of ef- 
fective pinning energy U* — T / S is convenient to extract 
this value. 

We can define inverse power-law form of flux activation 
energy U(J) as 



U(J) 



U 



[(Jco/JT i] 



(2) 
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Combining this with U — Tln(t/i e ff) extracted from Ar- 
rhenius relation, we can deduce the so-called "interpola- 
tion formula": 



J(T,t) 



■I 



cO 



[1 + (JJT/U ) Ht/tes) 



(3) 



where J c q is temperature dependent critical current den- 
sity in the absence of flux creep. From Eqs. (|2j and ([3]), 



U* = U + nTHt/tes) = U (J c0 /J) A 



(4) 



is derived. Thus the slope in double logarithmic plot of 
U* vs. 1/ J gives the value of fi, shown in the main panels 
of Fig. O In this way, we evaluate (i = 1.09 and 0.82 for 
pristine and H + -irradiated samples, respectively. Note 
that fj, ~ 1 in pristine crystal is often reported in Y-Ba- 
Cu-C»22 and IBSs 23 . Contrary to the above prediction 
of /i > 0, negative slope is observed at small J. This 
negative slope is often denoted as p in plastic creep the- 
ory with p = —0.5, and confirmed experimentally. 2 ^ Our 
evaluation of p — —0.61 in pristine sample is very similar 
to this value. 

To determine actual flux activation energy, we em- 
ploy extended Maley's method.— Since temperature de- 
pendence of U is not considered in the original Maley's 
method^ it is impossible to scale U in a wide range of 
J even if glassy exponent is unique. In order to solve 
this problem, appropriate temperature dependences of 
Uq and J C Q are assumed as follows. 



\2ln 



U (T) = U 00 [1-(T/T C ) 
J c0 (T) = J c0 o[l-(T/r c ) 2 ]". 



(5) 
(6) 



In order to simplify the problem, we choose the same ex- 
ponents in Eq. ([3]) and Eq. ©. Here exponent n is set to 
3/2, as in the case of Ref. WlM, while (1-T/T C ) 3 / 2 is se- 
lected in Bai_ x K :c Fe2As2i2 3 . U = -T\n[dM(t)/dt] + CT 
and C = \n(Buja/2iTr) is assumed as a constant, where 
B is the magnetic induction, uj is the attempt frequency 
for vortex hopping, a is the hopping distance, and r is 
the sample radius. We select C — 18 and 20 for pris- 
tine and H + -irradiated samples, respectively. Figures 
|3] show current density dependence of U in (a)pristine 
and (b)H + -irradiated Ba(Feo.93Coo.o7)2As2, respectively, 
constructed by extended Maley's method. Solid lines 
indicate power-law fitting to large J region where the 
slope in Fig. [5] is positive. Note that deviation of the 
data from the fitting in the small J region is reason- 
able since creep is plastic there. The obtained glassy 
exponents are (i = 1.01 and 1.24 for pristine and H + - 
irradiated samples, respectively. For pristine sample, this 
value is nearly the same as that obtained in Fig. [2ja), 
fi = 1.09. On the other hand, the change of fi by H + 
irradiation has an opposite trend. Namely, the value of 
[i decreases in Fig. [^b) , while grows in Fig. [SJb) after 
H + -irradiation. This is because the vortex system in H + - 
irradiated sample crossovers from elastic to plastic creep 
more gradually, as we can see in the main panel and in- 
set of Fig. [UJb). Hence we may underestimate /i and 
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FIG. 2. Inverse current-density dependence of effective 
pinning energy U* in (a)pristine and (b)H + irradiated 
Ba(Feo.93Coo.o7)2As2. Inset: Temperature dependence of 
normalized relaxation rate S. Solid line indicates fitting by 
Eq. Q. 



TABLE I. Parameters obtained from extended Maley's 
method and inverse of plateau value in S(T). 



Sample 

Pristine 
H + irradiated 



JcO 

(MA/cm 2 ) 
0.85 
2.90 



U 00 
(K) 
41.2 
93.1 



H nln(t/t eS ) 



1.01 
1.24 



35 
43 



overestimate p with the scheme of Fig. [21 and it is more 
reliable to estimate it from the extended Maley's method 
of Fig. |31 which uses more data points. For this reason, 
we conclude that p is slightly increased by H + irradia- 
tion. Additionally, slight increase of p is consistent with 
the regime of measurement, where it is closer to small 
bundle regime with p = 3/2, as we discussed in Fig. [T] 
Other resultant parameters are summarized in Table |U 
With these Uqo, temperature dependence of S is fitted 
by Eq. (p} with a single free parameter of plateau value 
gsat = i/n\ n (t/t cS ) as shown in the inset of Fig. H The 
inverse of this value is also shown in Table |U 

Using parameters obtained above, we calculate J after 
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FIG. 3. Current density dependence of U in (a) pristine and 
(b) H + -irradiated Ba(Feo.93Coo.o7)2As2 constructed by ex- 
tended Maley's method. Solid line indicates power-law fitting 
in large J region. 
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FIG. 4. Temperature dependence of J 3 (•), model function 
of J c before and after creep with the parameters in Table 
U (solid line), ST C - and <5Z-pinning (dashed line) and J c re- 
constructed by GIS (+) in (a)pristine and (b)H + -irradiated 
Ba(Feo.93Coo.o7)2As2. 



creep from attempt function of (true) critical current den- 
sity Eq. ([6]), which is shown as lower solid line in Fig. |4] 
In both cases, J is reasonably reproduced, especially at 
high J region (i.e. at low temperature). This means 
that the present collective pinning/creep analysis is ap- 
propriate. To get more insight into pinning mechanism 
in IBSs, we also show a function of ST C - and r5Z-pinning 
in Fig. |4] These functions are written as J c (t) / J c {0) = 
(l-i 2 ) 7/6 (l + * 2 ) 5/6 and (1 -i 2 ) 5 / 2 (l + i 2 )- 1 / 2 , respec- 
tively. 1 ' From this figure, our model function of J c can be 
considered as a superposition of the two pinning mech- 
anisms. To discuss such a mechanism, generalized in- 
version scheme (GIS) is utilized ! 27 ' 28 Although in this 
scheme, we have to assume empirical temperature de- 
pendence of penetration depth A and coherence length 
£ as cx (1 - t^y 1 ' 2 and oc (1 + t 2 )^ 2 {l - t 2 )- 1 ' 2 , re- 
spectively, we can directly reconstruct true critical cur- 
rent density J c from J s and discuss pinning mechanism. 
When we assume hi(t/t e s) ~ 23 in the measurement 
with field sweeping, and choose parameters for three- 
dimensional single vortex pinning, J c is reconstructed 
as shown in Fig. 2J which are in reasonable agreement 
with the model function. Similar analyses of pinning 



mechanism using GIS in pristine Ba(Fei_ :E Co :E )2As2 have 
been attempted in Ref. H^. They also conclude that 
both 5T C - and r5Z-mechanisms are working in this sys- 
tem. Here, we want to compare our work with similar 
work by Haberkorn et al.. 5 In their work, temperature 
dependence of (measured) current density J(T) is used 
to discuss the pinning mechanism in pristine and proton- 
irradiated Ba(Fei_ a; Co a ;)2As2. However, identification of 
pinning mechanism using J(T) is only empirical and lacks 
firm physical background. So, although their conclusion 
and our conclusion on the pinning mechanism are similar, 
we believe that our identification of pinning mechanism 
is more appropriate. 

We can basically describe physical quantities for vortex 
system by means of collective creep theory both in the 
case of pristine and H + -irradiated Ba(Feo.93Con.o7)2As2. 
It is noteworthy that the effect of H + irradiation can 
be summarized as E/no enhancement, without replacing 
model function Uo(T)/Uqo and Jo(T)/Joo- Namely, we 
can conclude that the effect of H + irradiation is enhance- 
ment of the collective pinning force by increasing weak- 
point-pinning centers without a drastic change of pinning 



5 



mechanism. 

Finally, we comment on the absolute value of J in- 
stead of J(T)/J(0). J is determined by the sum of 
weak-collective-pinning contribution J wc p and strong- 
point-pinning contribution J s pp.2^ Instead, if we assign 
J to J spp in the absence of flux creep (so that we write 
J as J c here), we can estimate upper limit of strong- 
pinning center fraction in the crystal. In the strong pin- 
ning theory^, critical current density is written as J c ss 
0.14^ 7 [DF(T)] 3 / 2 J . Here, J = c^ /l2^ 2 UKb « 
depairing current density, 7 = H^/H^ is anisotropy pa- 
rameter, n, D are density and diameter of pinning cen- 
ters, respectively. Assuming D as several times of £ a ;,, 
we can simplify F(T) m ln[l + D 2 /8$, 2 {T)} w 1. Using 
pinning center volume v « D 3 /2, nu = (J c /Jo) 2 (v2 x 
0.14 7 )- 2 w 0.05% with £ afc - 34 A from H c2 (0) ~ 280 
kOe£ and X a b ~ 2000 Ai22 This value is similar to the 
value reported in Na-doped CaFe2As2^ 

In summary, we have studied the effect of proton ir- 



radiation up to 1.2 x 10 16 cm~ 2 in optimally Co-doped 
BaFe2As2 single crystals. Critical current density under 
self-field is enhanced by a factor of 2.5 at 2 K. Temper- 
ature dependence of critical current density and normal- 
ized flux relaxation rate is interpreted by collective creep 
theory. With Maley's method, a glassy exponent /Lt ~ 1 
and variation of barrier height for flux creep Uoo ~ 41 
K to 93 K are directly determined. To explain the value 
of J after the creep from the model function of J c , J c is 
concluded to be controlled by both ST C ~ and 51— pinning. 
This model function is consistent with the result of gen- 
eralized inversion scheme. Proton irradiation effect is 
concluded as doubling of barrier height in the absence of 
flux creep. 
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